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Abstract. The chemical composition of a protoplanetary disk is determined not only by in situ 

!_. . chemical processes during the disk phase, but also by the history of the gas and dust before 

it accreted from the natal envelope. In order to understand the disk's chemical composition at 

the time of planet formation, especially in the midplane, one has to go back in time and retrace 

the chemistry to the molecular cloud that collapsed to form the disk and the central star. Here 

i i' we present a new astrochemical model that aims to do just that. The model follows the core 

•^T collapse and disk formation in two dimensions, which turns out to be a critical upgrade over 

rh ■ older collapse models. We predict chemical stratification in the disk due to different physical 

conditions encountered along different streamlines. We argue that the disk-envelope accretion 

_^ shock does not play a significant role for the material in the disk at the end of the collapse phase. 

Finally, our model suggests that complex organic species are formed on the grain surfaces at 

temperatures of 20 to 40 K, rather than in the gas phase in the T > 100 K hot corino. 
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^. ■ 1. Introduction 

We live on a world rich in complex organic molecules, all of which ultimately derive from 
{-*. a cold pre-stellar core with a much simpler chemical composition. How did the chemistry 

jvq ■ in the solar system evolve from that original simple state to its current complex state? 
j^J This chapter focuses on one particular part of the entire evolutionary track: the transition 

r-^ ■ from a pre-stellar core to a protoplanetary disk. 

y—{ ' The chemical evolution in a collapsing envelope has been studied extensively with one- 

r—( • dimensional (ID) spherical models (e.g., Bergin & Langer 1997; Rodgers & Charnley 
2003; Lee et al. 2004; Aikawa et al. 2008), but such models cannot properly describe 
• ■ *h the circumstellar disk. We have recently constructed a new collapse model that describes 

the envelope and the disk in two dimensions (2D) as an axisymmetric structure (Visser 
et al. 2009a, 2011; Visser & Dullemond 2010). By taking the vertical extent of the disk 
into account, we find that accretion from the envelope onto the disk largely takes place 
near the cold outer edge, rather than near the warmer inner parts as suggested by earlier 
models. The lower temperatures in our scenario allow for a larger fraction of pristine 
material to enter the disk. The 2D models also show strong evidence of layered accretion 
in the disk: new material arriving from the envelope always accretes on top of the older 
material already in the disk (Fig. Q]). Another improvement in our model is that we 
calculate the dust temperatures throughout the disk and the envelope with a full radiative 
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Figure 1. Schematic view of layered accretion. The outer parts of the original cloud core (left) 
end up as the surface layers of the circumstellar disk (right). The white part of the core is 
roughly the part that disappears into the outflow. The arrows illustrate the flow of material 
onto the disk. If the vertical extent of the disk is not taken into account properly, as happened 
in some older models, such streamlines would extend further in and the radii at which material 
accretes would be underestimated. 

transfer code. The approximate temperature calculations used by some older models are 
not accurate enough for a reliable investigation of the chemistry. 

We provide here a simple summary of our model (Sect. [5]) and highlight some of the 
first results (Sect. \S§. In particular, we focus on the chemical history of H 2 and other 
simple molecules as material flows from the pre-stellar core to the planet-forming disk. 
The discussion (Sect. 0| addresses two open questions in star formation: how strong is 
the disk-envelope accretion shock, and do hot corinos really exist? 



2. Model description 

In order to follow the chemical evolution from a pre-stellar core to a circumstellar disk, 
we need a physical framework to provide densities, temperatures and other quantities that 
may affect the chemistry. We choose to use a fast semi-analytical model for the density 
and velocity profiles and a full continuum radiative transfer code for the dust temperature 
and the radiation field. The dust temperature is particularly important for the gas-grain 
interactions. Hydrodynamical codes often use an approximation of the dust temperature 
and are therefore less suited to study the chemical evolution. 

The full model procedure is explained in detail in Visser et al. (2009a), Visser & 
Dullcmond (2010) and Visser et al. (2011). We summarise the model in this section by 
way of the step- wise description from Fig. [2] 

2.1. Initial conditions 

The physical and chemical evolution of the collapsing envelope depends strongly on the 
initial conditions. The three most prominent parameters in our model are the mass of 
the cloud core before the onset of collapse (Mo), the isothermal sound speed (c s ), and 
the core's solid-body rotation rate (Jlo). A higher rotation rate leads to a larger and 
more massive disk, as do a lower sound speed and a higher core mass (Visser et al. 
2009a). Lowering the sound speed or increasing the core mass both tend to produce 
a colder disk, with larger fractions of ices relative to gas-phase species. Increasing the 
solid-body rotation rate also leads to a colder disk. The core mass and the sound speed 
together determine how long it takes for the entire core to accrete onto the star and disk: 
t acc = M /M a M G/cl (Shu 1977). 
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2.2. Density and velocity 

We follow the inside-out solution of Shu (1977) for a gravitationally collapsing spherical 
envelope. However, because of the core's rotation, spherical symmetry is lost at the onset 
of collapse. The effects of rotation can be treated as a perturbation of the pure spherical 
solution (Cassen & Moosman 1981; Terebey et al. 1984). In the perturbed solution, 
matter no longer streams in radially, but it is deflected towards the gravitational midplane 
and builds up a circumstellar disk. After the disk is first formed, it evolves by ongoing 
accretion from the envelope and by viscous spreading to conserve angular momentum 
(Shakura & Sunyaev 1973; Lynden-Bell & Pringle 1974; Dullemond et al. 2006). 

The size of the disk is initially set by the centrifugal radius (R c w (l/16)Cg£ 3 ilo, with 
t the time since the onset of collapse), because the streamlines from the envelope only 
hit the midplane inside of that point. However, viscous spreading can push the outer 
edge of the disk to several times R c (Visser et al. 2009a; van Weeren et al. 2009). We 
adopt a Gaussian profile for the vertical density structure of the disk, from which we also 
derive the vertical velocity profile. The boundary between the disk and the envelope is 
conveniently defined as the surface where the thermal pressure from the disk equals the 
ram pressure from the infalling envelope (Visser & Dullemond 2010). 

The analytical collapse solution of Terebey et al. (1984) does not include a bipolar 
outflow, so we have to put one in as an ad-hoc addition. Based on observations and 
theoretical predictions (Velusamy & Langer 1998; Canto et al. 2008), we adopt curved 
cavity walls whose outflow angle grows with time; specifically, at a height z above the 
midplane, the size i? C av of the cavity is proportional to t 2 z 2 / 3 (Visser et al. 2011). 

2.3. Temperature and radiation field 
The chemical evolution depends strongly on the temperature, so we employ a full radia- 
tive transfer method instead of an analytical approximation for step 3 of our model. The 
only heating source taken into account is the central star, whose luminosity and effective 
temperature evolve according to Adams & Shu (1986) and Young & Evans (2005). Tak- 
ing the density profiles from step 2, we compute the 2D dust temperature profiles at a 
range of time steps with RADMC (Dullemond & Dominik 2004). As additional output, 
RADMC also provides the intensity of the radiation field as function of wavelength at 
each point in the disk and remnant envelope, as required for the photon-driven reactions 
in our chemical network. 

We set the gas temperature equal to the dust temperature throughout the disk and 
the envelope, even though this is a poor assumption in the surface of the disk and the 
inner parts of the envelope (Kamp & Dullemond 2004; Jonkheid et al. 2004). Planned 
updates to the model include a separate treatment of the dust and gas temperatures. 
The effect of the accretion shock on the temperature profiles is discussed in Sect. |U 

2.4. Infall trajectories 

We follow the chemical evolution as material streams from the natal core towards the 
circumstellar disk, so we have to solve the chemistry in a Lagrangian frame. The semi- 
analytical nature of the density and velocity solutions offers two ways to approach the 
problem. The first way, followed in Visser et al. (2009a) , consists of populating the core 
with a large number of points (at least 10000) and calculating the infall streamlines in 
the forward direction. Alternatively, as done in Visser et al. (2011), we can select points 
of interest in the disk at the end of the collapse phase and calculate the trajectories 
backwards in time to their point of origin in the core. In both cases, we end up with a set 
of trajectories (Cartesian or polar coordinates as function of time) with corresponding 
physical conditions (density, temperature, UV field and extinction as function of position 
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Figure 2. Step- wise summary of our 2D axisymmetric collapse model. Steps 2 and 4 are semi- 
analytical, while steps 3 and 5 consist of detailed numerical simulations. (Figure adapted from 
Visser et al. (2011).) 

and time). These are used as input for the chemistry code. Each trajectory is treated 
in isolation, i.e., the material streaming along one trajectory is assumed not to interact 
with material streaming along other trajectories. 

2.5. Chemistry 

The basis of our chemical network is the UMIST06 database (Woodall et al. 2007) as 
modified by Bruderer et al. (2009), except that X-ray chemistry is not included. Full 
details are provided in Visser et al. (2011). The network includes freeze-out and evap- 
oration (both thermal and photon-driven), but we limit the grain-surface chemistry to 
a handful of simple hydrogenation reactions: atomic H, C, N, O and S can react with 
incoming H to form H2, CH4, NH3, H2O and H2S (Black & van Dishoeck 1987; Bergin 
& Langer 1997; Hollcnbach et al. 2009). Other potentially important grain-surface reac- 
tions, such as the formation of H2O out of O2 (Ioppolo et al. 2008) or the formation of 
H2CO and CH3OH out of CO (Fuchs et al. 2009), have not yet been incorporated. The 
gas-phase network includes photodissociation and photoionisation reactions (accounting 
specifically for self-shielding of CO; Visser et al. 2009b), which dominate the chemistry 
close to the outflow wall and in the surface layers of the disk. 

In order to set the chemical composition at the onset of collapse (t = 0), we evolve the 
initially atomic gas for a period of 1 Myr at nn = 8 x 10 4 cm~ 3 and T gas = Td us t = 10 
K. The resulting solid and gas-phase abundances are consistent with those observed in 
pre-stellar cores (e.g., di Francesco et al. 2007). They form the initial conditions for the 
chemical evolution along each infall trajectory. 



3. Results 

3.1. Key chemical trends 

This section highlights a few key results from the standard model of Visser et al. (2009a, 
2011), which has initial conditions Mo — 1.0 M Q , c s = 0.26 km s _1 and £!o 
The chemistry is run for the full duration of the envelope accretion phase (i ac 
yr), following the 1 Myr pre-collapse phase. 
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Figure 3. Schematic view of the chemical evolution along one infall trajectory (solid black 
curve). The solid and dashed grey lines denote the surface of the disk and the outflow wall, both 
at t — t acc . The abundances of some key species (black bars: gas; grey bars: ice) are plotted at 
four points along the trajectory (A: top panel; B: right panel; C: left panel; D: bottom panel). 
The key processes governing the overall chemistry at each point are listed in the bottom right. 
(Figure adapted from Visser et al. (2011).) 



We first focus on the chemical evolution along one particular trajectory, starting near 
the edge of the original core (i?; = 5050 AU, z\ = 4420 AU) and ending at t = t acc in 
the inner part of the disk, about 0.2 AU below the surface (R{ — 6.3 AU, Z{ = 2.4 AU). 
The chemistry along this and other trajectories is governed by changes in the physical 
conditions. Figure |3] shows how the solid and gas-phase abundances of eight important 
species change along the streamline. We identify four chemical events: 

• At point A, the temperature reaches a value of about 20 K and simple ices like CO, 
N2 and O2 begin to evaporate. 

• The temperature increases further to 30-40 K at point B, driving CH4 and NO off 
the grains. 

• Moving in further, the dust is heated to the evaporation temperature of H2O (100 
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Figure 4. Examples of qualitatively different infall trajectories terminating at nearly the same 
point at the end of the accretion phase: R — 2.9 AU (top) and 3.7 AU (bottom). The black and 
grey portions of the trajectories indicate where most of the water is frozen out or in the gas 
phase. The points where the trajectories reverse direction are marked with the time in years. 

K). This close to the outflow cavity wall, the gas is exposed to the stellar UV field, 
dissociating the evaporating H2O and other previously evaporated species. Our model 
does not include an excess UV field for the young star, so there are not enough high- 
energy photons to dissociate CO. 

• Finally, at point D, the high density (10 10 -10 11 cm~ 3 ) allows some of the dissociated 
species to reform in ion-neutral and neutral-neutral reactions. 

3.2. Chemical stratification 

Different streamlines exhibit different physical conditions. In order to fully understand 
the chemical evolution from every point in the pre-stellar core to every point in the 
circumstellar disk, we have to investigate more than just the single trajectory from Fig. 
G3 In this regard, it is interesting to look at where material from different parts of the 
core ends up in the disk. 

The hydrodynamical simulations of Brinch et al. (2008) and our semi- analytical model 
(Visser et al. 2009a) both show the envelope to accrete onto the disk in a layered fashion 
(Fig. [1]). Material from near the centre of the pre-stellar core is the first to become part 
of the disk, and material from larger radii continuously piles up on top. The thickness 
of the disk is small relative to its radial size (typical scale heights are 0.1-0.2i?), but the 
outer edge is high enough that a substantial fraction of the infalling material hits there 
rather than on the surface (Fig. [T] see also Fig. 3 in Visser & Dullemond 2010). The 
infalling material has to compete with the outwardly expanding disk, resulting in some 
material being pushed up and over the outer edge in the hydrodynamical simulations 
(see also van Weeren et al. 2009). 

Individual streamlines can exhibit large qualitative differences despite terminating at 
nearly the same point in the disk. Figure [4] shows two such streamlines, with termini at 
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Figure 5. Schematic view of the history of H2O gas and ice towards different parts of the disk. 
The typical gas-phase abundance is indicated for each zone (zones 4 and 5 both have 10 -4 ). The 
position of the snowline (R = 5 AU, between zones 5 and 6) is marked with an asterisk. Note 
the disproportionality of the R and 2 axes. (Figure adapted from Visser et al. (2011).) 



2.9 and 3.7 AU from the star. The streamline in the bottom panel originates at r = 1100 
AU from the centre of the pre-stellar core. It reaches the disk at a fairly early time 
(4.2 x 10 4 yr after the onset of collapse), when the disk is only a few AU large. The 
parcel of gas following this streamline gets swept up in the outward viscous expansion 
until, at t = 1.7 x 10 5 yr, enough mass has accreted at larger radii that it gets pushed 
back in again. The streamline in the bottom panel originates further out in the natal 
core, at r = 3500 AU, and therefore accretes at a later time, t = 1.3 x 10 5 yr, when the 
disk has grown to about 20 AU. It experiences a brief period of outward motion before 
moving in further. 

In terms of the chemical evolution of the material flowing along these two streamlines, 
the key difference lies in the conditions encountered before t — 10 5 yr. The bottom 
streamline passes much closer to the star than does the other one (1.6 versus 16 AU); it 
is heated to more than 500 K, whereas the top parcel initially does not get warmer than 
120 K. The latter is enough to drive H2O and other non-refractory species off the grains 
(as indicated by the grey portion of the streamline up to 1.3 x 10 5 yr), but not enough 
to power the kind of high-temperature gas-phase chemistry that is possible at 500 K. 

The different types of streamlines give rise to different chemical zones in the disk. This 
effect of "chemical stratification" is visualised in Fig.[5j which traces the history of H2O 
gas and ice towards various parts of the disk. (Similar figures for CH4 and NH3 appear in 
Visser et al. (2011).) The streamlines terminating in zone 1 never encounter temperatures 
above 100 K, so H2O always remains frozen on the grains. Zone 2 contains the terminus 
of the streamline from Fig. [3] some H2O ends up here, but most of it is dissociated by 
the stellar radiation field. The same happens for material ending up in zone 3, except 
that no H2O reforms at all after being destroyed. Parcels ending up in zone 4 do not 
encounter a strong radiation field; H2O in these parcels stays on the grains until passing 
the snowline at ~5 AU. Finally, zones 5 and 6 consist largely of material that accreted 
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at early times, following trajectories like the one in the bottom panel of Fig. 2) H2O in 
these parcels undergoes a freeze-thaw cycle: it evaporates during the initial inwards part 
of the trajectory and freezes out again during the outwards leg. Parcels ending up in 
zone 5 undergo another evaporation event as they pass the snowline; material in zone 6 
lies beyond the snowline, so H 2 remains frozen the rest of the time. 

The chemical stratification comes with one major caveat: mixing. In our model, there 
is no exchange of material between adjacent parcels. In reality, however, turbulent mixing 
is likely to change the sizes and locations of the chemical zones, and to make the borders 
between the zones more diffuse (e.g., Semenov et al. 2006). Spectroscopic observations at 
5-10 AU resolution, as possible with ALMA, are required to determine to what extent 
chemical stratification occurs in actual disks. 



4. Discussion 

4.1. How strong is the disk-envelope accretion shock? 

When the supersonically infalling material from the envelope hits the nearly static disk, 
it passes through a J-type shock. However, this shock may not play as large a role as 
commonly believed. Neufeld & Hollenbach (1994) showed that the strength of the shock 
decreases with decreasing pre-shock density and decreasing shock velocity. For a given 
radius, the density and velocity decrease as the collapse proceeds; for a given point in 
time, both decrease with radius (Cassen & Moosman 1981; Terebey et al. 1984). 

During the second half of the collapse phase (from t = 0.5i acc to t acc ), at least 80% of 
the infalling material accretes onto the outer half of the disk. (The exact number depends 
on the initial conditions.) At the end of the collapse phase, the bulk of the disk consists 
of material that accreted at fairly late stages and large radii. Most of the material that 
accreted at early times or small radii gets pushed into the star. Hence, we expect that 
the bulk of the gas and the dust in a T Tauri or Herbig Ae/Be disk did not pass through 
a strong accretion shock. Quantitatively, we showed in Visser et al. (2009a) that the 
shock has a smaller effect on the temperature history of the disk material than does 
the radiation from the young star. In terms of the chemical evolution, this means the 
desorption of H2O and organic species is likely unaffected by the accretion shock. Species 
with lower binding energies, however, may be more sensitive. 

4.2. Do hot corinos exist? 

Observations of complex organic species like methyl formate (HCOOCH3) and dimethyl 
ether (CH3OCH3) in low-mass protostars have been linked to the high-temperature 
(> 100 K) inner envelope, also known as the hot core or hot corino (Ceccarelli 2004). 
One possible formation pathway is powered by the evaporation of H2O and other non- 
refractory species at 100 K, followed by a chain of high-temperature gas-phase reactions 
to transform simple carbon-bearing molecules into second-generation complex organics 
(Herbst & van Dishoeck 2009; see also Charnley et al. 1992). Alternatively, the organics 
may be formed abundantly on the grain surfaces through mild photochemistry at 20-40 
K (Garrod & Herbst 2006; Garrod et al. 2008). Which scenario is more likely? 

The first formation route requires the gas to spend at least several 10 3 yr at temper- 
atures above 100 K (Charnley et al. 1992). In our model, however, the material in the 
inner envelope is essentially in freefall towards the star or the inner disk, and it only 
spends a few 100 yr at the required temperatures (Visser et al. 2009a; see also Schoier et 
al. 2002). The grain-surface pathway takes several 10 4 yr at 20-40 K (Garrod & Herbst 
2006). Such temperatures are attainable both in the envelope before the material reaches 
the hot inner part and inside the disk (e.g., Fig. 12 in Visser et al. 2009a). Once formed 
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on the grains, these first-generation organics can be returned to the gas phase by thermal 
evaporation in the hot corino, the hot inner disk, or the surface layers of the outer disk. 
They can also be liberated by non-thermal processes such as photodesorption in the disk 
and shocks in the outflow walls. In any scenario, the mass of the hot corino is quite lim- 
ited: the amount of material wedged between the disk and the outflow cavity is at most 
a few per cent of the amount of material inside the disk. Hence, even if the hot corino is 
rich in organics, there still would not be enough of them to explain the observations. 

Based on the timescales and the mass, we conclude that the observed second-generation 
complex organics are not associated with the hot core or hot corino. The circumstellar 
disk and outflow walls offer a more likely explanation: they have a larger mass reservoir 
and are dynamically stable enough for the observed organics to form abundantly. 



5. Conclusions 

We have constructed a two-dimensional axisymmetric model to simulate the collapse 
of a pre-stellar core into a young star and its surrounding disk. We use the model to 
study the chemical evolution of the gas and the dust streaming from the envelope to 
the disk. Comparing the streamlines to what happens in one-dimensional models, we 
conclude that it is imperative to study core collapse and disk formation in 2D. Our 
model predicts chemical stratification in the disk due to different physical conditions 
(density, temperature, radiation field) encountered along different streamlines. We argue 
that the disk-envelope accretion shock does not play a significant role for the material 
in the disk at the end of the collapse phase. Finally, our model suggests that complex 
organic species are formed on the grain surfaces at temperatures of 20 to 40 K, rather 
than in the gas phase in the T > 100 K hot corino. 
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Discussion 

Yamamoto: How do you assume the initial conditions in your model? We now know 
that various complex molecules exist in protostellar cores. How are they brought into the 
inner disk? 

Visser: We have a static pre-collapse phase of 10 6 yr to produce typical pre-stellar 
core abundances. Since we do not have any grain-surface chemistry other than simple 
hydrogenation of O, C, N and S, complex organics are not formed abundantly. However, 
based on the models by Garrod & Herbst (2006) and Garrod et al. (2008), we predict 
that complex organics are formed abundantly inside the disk. We have not considered 
warm carbon-chain chemistry. 

Johnstone: How important are accretion shocks for the chemical evolution of infalling 
gas? 

Visser: The accretion shock is not important for the bulk of the material in the disk 
at the end of the accretion phase, because that material accreted at late times and 
large radii, and therefore did not encounter a strong shock (Neufeld & Hollenbach 1994). 
Material that accretes earlier encounters a stronger shock, but this material disappears 
into the star while the rest of the disk is formed. We refer to Sect. 14. ll for more details. 

AlKAWA: What is the cause of the vortices that appear in the animation? Do they also 
appear in your 2D analytical model? 

Visser: The vortices are due to material piling up at the outer edge of the disk and 
being pushed up and over the outer edge in the 2D simulation. They do not appear in 
the semi-analytical model. 

Irvine: Evidence suggests that most stars, including the Sun, did not form in isolation 
(the Shu model). Have you considered how this might affect your results? 

Visser: No, we have not considered this in any detail. We expect that many of the 
results hold for circumbinary envelopes and disks, as long as the interbinary distance is 
smaller than 10 AU or so. If the protostar in question forms in the vicinity of an O or 
B star, the strong external radiation field would result in a larger fraction of atomic and 
ionised species at the outer edge of the envelope and the disk. 



